(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     28453,        711]
NotebookOptionsPosition[     26761,        654]
NotebookOutlinePosition[     27115,        670]
CellTagsIndexPosition[     27072,        667]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[StyleBox["At T = 0.12 T_F :",
 FontColor->RGBColor[1, 0, 0]]], "Text",
 CellChangeTimes->{{3.804734401936799*^9, 3.8047344186024942`*^9}, {
  3.819723761863247*^9, 3.819723771136421*^9}, {3.819723824275384*^9, 
  3.819723846937909*^9}, {3.819723898656394*^9, 3.819723913256065*^9}, {
  3.820461950722602*^9, 3.820461973085478*^9}, {3.820636715442717*^9, 
  3.8206367526214523`*^9}, {3.821740366958331*^9, 3.821740368075968*^9}, {
  3.821740399955288*^9, 3.821740402686111*^9}, {3.821740636393045*^9, 
  3.82174065714289*^9}, {3.821901958417396*^9, 3.8219019653092117`*^9}, {
  3.822441473362677*^9, 3.82244149834795*^9}}],

Cell[TextData[StyleBox["With 10 * 19400 atoms: Recoil energy 518 nK & Fermi \
energy 455 nK & T = 54.6 nK (for Fig. 4)",
 FontColor->RGBColor[1, 0, 0]]], "Text",
 CellChangeTimes->{{3.804734401936799*^9, 3.8047344186024942`*^9}, {
  3.819723761863247*^9, 3.819723771136421*^9}, {3.819723824275384*^9, 
  3.819723846937909*^9}, {3.819723898656394*^9, 3.819723913256065*^9}, {
  3.820461950722602*^9, 3.820461973085478*^9}, {3.820636715442717*^9, 
  3.8206367526214523`*^9}, {3.821740366958331*^9, 3.821740368075968*^9}, {
  3.821740399955288*^9, 3.821740402686111*^9}, {3.821740636393045*^9, 
  3.82174065714289*^9}, {3.821740694573724*^9, 3.821740710730578*^9}, {
  3.821741627887662*^9, 3.821741628945161*^9}, {3.821901784306122*^9, 
  3.82190178700977*^9}, {3.821901916652829*^9, 3.821901953841399*^9}, {
  3.822441505892055*^9, 3.822441506241353*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"mm", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"10", "^", "18"}], ")"}], " ", "87", " ", 
    RowBox[{"166", "/", "100"}], " ", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "27"}], ")"}]}]}]}], ";", 
  RowBox[{"wx", "=", 
   RowBox[{"2", " ", "Pi", " ", 
    RowBox[{"120", "/", "1.0"}]}]}], ";", 
  RowBox[{"wy", "=", 
   RowBox[{"2", " ", "Pi", " ", "120", "*", "1.0"}]}], ";", 
  RowBox[{"wz", "=", 
   RowBox[{"2", " ", "Pi", " ", "506"}]}], ";", 
  RowBox[{"who", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"8", " ", 
      RowBox[{"Pi", "^", "3"}], " ", "506", " ", "120", " ", "120"}], ")"}], 
    "^", 
    RowBox[{"(", 
     RowBox[{"1", "/", "3"}], ")"}]}]}], ";", 
  RowBox[{"kB", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"10", "^", "6"}], " ", 
      RowBox[{"10", "^", "6"}], " ", 
      RowBox[{"10", "^", "18"}]}], ")"}], " ", "1.38", " ", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "23"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJwd0n0s1HEAx/G7kxxJRlOhojTSEatoiH1rZ9clVKR56FSiR2LjxIrqehij
B01CnopQOCm1prqk6wilVohapUf3TXmmVev3+f7x2eu/z19v6+0xG3cKeDye
mBvcmCL33vBIS6Rh48zVhc+lUNFqFgSr9MQymCa4tRvWK17KYVG7VyIU7i65
HsjpYi1RwqIbe+4e5NSE1TCzy4Rq2KDn/AQGzAuwbXisJd119Uyjt7PudWu1
RBT1/T7MTR59AlP8AluhcsCsA5YFVXfBX6WidzCwP5rpoWrJuWBMSaT2aDGc
ccdxYb4VJR5LLW3gAkHBynI7Sl5EK5lj6+fpS9wouS2WM/3tXU2gm/sdpoWg
xQzyzC1mwVp/y0MjWyhZp1jL9Ciry4eOPfxLUN0aUgpPSBOYzfM7U2tDKNkU
8i0Tqp2Wuk/uoOShuMMT5jQ+9YU3ht4yZ9f8/OUaQ4m9dOcITMwL7mzhNF7n
0wMN5o7JphygxEpjEw7DU7z2w66I5dHweaZTKfRp28X8klwxTZcz2EjDPJN5
LtaQM6P8LjNH/6IC2p5XMT8OT2TAzy4LM6HiSNVleJL2MWdW52quch5JamR6
142MDnCu6RSNw+Pt4qPOsZT03ithnspfXQz1tx1j8iITKuG17AfMtMOrbkJP
mZypXkVVUGX64yF8YH6qFT4LznoGe2Sejss5Iz7HOcOATmH3kjhK3F+ZvIGR
4+e+wa21r5ld8RWDIs6YSoMhGNBs/gemZkz7C4WqDr4D5+mmFTpQ/29EYVo8
JXHWViUwkq+rhMdVEmZqzMB9aNTb1Ajtjkk0MFaS3gJfDzZ/gFTw7xMUqq/y
0/FTMMRU9Yvs4ItFB5mXTV+GQd1FS8Lh3u2DO2BO6PQI6K96vxdW1Ibvh1Tx
tRLuMfGrgg1lhT6jnNH1832h1VRlPNTRJiVA24kVnjYJXG99fC/ooGe9eBnn
rX1nHOBAk7wIRg26FsPQdKHUhdNgsoG5ayyqd5hzztnEdzDL0LYP2m+WMZ3s
rnzQyaAkaWv5R+gRWjHelkfJm2H+b/gfbFE3zQ==
  "]],

Cell["position in 10^-6 m (um), mass in 10^-18 kg (femtogram)", "Text",
 CellChangeTimes->{{3.8197216308824797`*^9, 3.8197216454928207`*^9}, {
  3.819725875776466*^9, 3.8197258843290367`*^9}, {3.819940595715087*^9, 
  3.8199405997537394`*^9}, {3.819940744608406*^9, 3.819940745946722*^9}, {
  3.820068611129483*^9, 3.820068615020266*^9}, {3.820068754146092*^9, 
  3.8200687611136723`*^9}, {3.8200688973296537`*^9, 3.820068907017367*^9}, {
  3.8397827741557283`*^9, 3.839782774682128*^9}, {3.8397828073798637`*^9, 
  3.8397828078341417`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"fugacity", "[", "tovertf_", "]"}], ":=", 
   RowBox[{"z", "/.", 
    RowBox[{"FindRoot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"PolyLog", "[", 
        RowBox[{"3", ",", 
         RowBox[{"-", "z"}]}], "]"}], "+", 
       RowBox[{"1", "/", 
        RowBox[{"(", 
         RowBox[{"6", " ", 
          RowBox[{
           RowBox[{"(", "tovertf", ")"}], "^", "3"}]}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "1000"}], "}"}]}], "]"}]}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fugacity", "[", "0.12", "]"}]], "Input",
 CellChangeTimes->{{3.819721573451778*^9, 3.8197215817604637`*^9}, 
   3.8200381415337*^9, 3.820066842720181*^9, {3.820066910479497*^9, 
   3.820066910868079*^9}, {3.820067319725153*^9, 3.820067320161503*^9}, {
   3.8222695256946983`*^9, 3.822269526073217*^9}, {3.8222729325233192`*^9, 
   3.822272933347632*^9}, {3.8320290365463667`*^9, 3.8320290368822737`*^9}}],

Cell[BoxData["2804.127417323298`"], "Output",
 CellChangeTimes->{
  3.81972158322726*^9, 3.819940922923462*^9, 3.819941217438712*^9, {
   3.820038145804867*^9, 3.820038152373678*^9}, 3.820066748240443*^9, 
   3.820066850825206*^9, 3.820066911768078*^9, 3.820067320769437*^9, 
   3.8200681473698797`*^9, 3.820073607185788*^9, 3.8204621512887783`*^9, 
   3.822269378167272*^9, 3.8222695271265497`*^9, 3.8222729346262712`*^9, 
   3.8320290382061768`*^9, 3.8320313024189568`*^9, 3.8397714576634827`*^9, 
   3.839783189321555*^9}]
}, Open  ]],

Cell["occupations :", "Text",
 CellChangeTimes->{{3.8397948338497887`*^9, 3.839794855285736*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"fuga", "=", "2804"}], ";", 
  RowBox[{"pxtrans", "=", "0.0006"}], ";", 
  RowBox[{"beta", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"kB", " ", "54.6", " ", 
      RowBox[{"10", "^", 
       RowBox[{"(", 
        RowBox[{"-", "9"}], ")"}]}]}], ")"}]}]}], ";", 
  RowBox[{
   RowBox[{"nfdi", "[", 
    RowBox[{
    "px_", ",", "py_", ",", "pz_", ",", "qx_", ",", "qy_", ",", "qz_"}], 
    "]"}], "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{
       RowBox[{"Exp", "[", 
        RowBox[{"beta", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wx", "^", "2"}], " ", 
            RowBox[{"qx", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wy", "^", "2"}], " ", 
            RowBox[{"qy", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wz", "^", "2"}], " ", 
            RowBox[{"qz", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{"px", "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}], "+", 
           RowBox[{
            RowBox[{"py", "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}], "+", 
           RowBox[{
            RowBox[{"pz", "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}]}], ")"}]}], "]"}], "/", 
       "fuga"}]}], ")"}]}]}], ";", 
  RowBox[{
   RowBox[{"nfdf", "[", 
    RowBox[{
    "px_", ",", "py_", ",", "pz_", ",", "qx_", ",", "qy_", ",", "qz_"}], 
    "]"}], "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{
       RowBox[{"Exp", "[", 
        RowBox[{"beta", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wx", "^", "2"}], " ", 
            RowBox[{"qx", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wy", "^", "2"}], " ", 
            RowBox[{"qy", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{"mm", "/", "2"}], " ", 
            RowBox[{"wz", "^", "2"}], " ", 
            RowBox[{"qz", "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"px", "+", "pxtrans"}], ")"}], "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}], "+", 
           RowBox[{
            RowBox[{"py", "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}], "+", 
           RowBox[{
            RowBox[{"pz", "^", "2"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", " ", "mm"}], ")"}]}]}], ")"}]}], "]"}], "/", 
       "fuga"}]}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.819721668898777*^9, 3.8197217061866817`*^9}, {
   3.81972175086093*^9, 3.819721857025168*^9}, {3.81972189373664*^9, 
   3.819721898184353*^9}, {3.8197242054309187`*^9, 3.819724207207218*^9}, {
   3.819724238740638*^9, 3.8197242808207493`*^9}, {3.819724311646991*^9, 
   3.819724327410256*^9}, {3.819724510852416*^9, 3.819724583404632*^9}, {
   3.819724654885054*^9, 3.819724689629346*^9}, {3.819724818523426*^9, 
   3.8197248799330587`*^9}, 3.819725348272904*^9, {3.819726141124792*^9, 
   3.819726145074005*^9}, {3.8197271925912046`*^9, 3.819727194520245*^9}, {
   3.819727268331959*^9, 3.819727274775401*^9}, {3.819727335365452*^9, 
   3.819727362360982*^9}, {3.820067350732499*^9, 3.82006735113525*^9}, {
   3.8200707406790648`*^9, 3.820070748067524*^9}, {3.820073612882368*^9, 
   3.820073633158723*^9}, {3.8204618784743967`*^9, 3.82046187930412*^9}, {
   3.822269534409564*^9, 3.822269559567038*^9}, {3.822272922635233*^9, 
   3.822272944122801*^9}, {3.832029069443962*^9, 3.8320290720486507`*^9}, {
   3.832029161448626*^9, 3.832029164473154*^9}, {3.832031319337771*^9, 
   3.8320313218034773`*^9}}],

Cell["\<\
spatially resolved suppression factor vs. radial distance qr :\
\>", "Text",
 CellChangeTimes->{{3.839790227987015*^9, 3.83979025252203*^9}, {
  3.839794964561625*^9, 3.839794971337887*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"nint01", "[", "qr_", "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"nfdi", "[", 
      RowBox[{"px", ",", "py", ",", "pz", ",", "qr", ",", "0", ",", "qz"}], 
      "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"px", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"py", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"pz", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"qz", ",", 
       RowBox[{"-", "10"}], ",", "10"}], "}"}]}], "]"}]}], ";", 
  RowBox[{
   RowBox[{"nint02", "[", "qr_", "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"nfdi", "[", 
       RowBox[{"px", ",", "py", ",", "pz", ",", "qr", ",", "0", ",", "qz"}], 
       "]"}], " ", 
      RowBox[{"nfdf", "[", 
       RowBox[{"px", ",", "py", ",", "pz", ",", "qr", ",", "0", ",", "qz"}], 
       "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"px", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"py", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"pz", ",", 
       RowBox[{"-", "0.005"}], ",", "0.005"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"qz", ",", 
       RowBox[{"-", "10"}], ",", "10"}], "}"}]}], "]"}]}], ";", 
  RowBox[{
   RowBox[{"sfactor", "[", "qr_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"nint01", "[", "qr", "]"}], "-", 
      RowBox[{"nint02", "[", "qr", "]"}]}], ")"}], "/", 
    RowBox[{"nint01", "[", "qr", "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.8197469728230047`*^9, 3.819746988377232*^9}, {
   3.8197471576181383`*^9, 3.8197471990477324`*^9}, 3.8197472714304028`*^9, 
   3.8197474338222523`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sfactorlist", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"qr", ",", 
      RowBox[{"sfactor", "[", "qr", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"qr", ",", "0", ",", "15", ",", "0.5"}], "}"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.8197495243497353`*^9, 3.819749571978118*^9}, {
  3.819749602308264*^9, 3.81974961258428*^9}, {3.819754706005807*^9, 
  3.819754724943501*^9}}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","slwcon",
   "\"Numerical integration converging too slowly; suspect one of the \
following: singularity, value of the integration is 0, highly oscillatory \
integrand, or WorkingPrecision too small.\"",2,6,1,21519474957897888619,
   "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783259326583*^9}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","eincr",
   "\"The global error of the strategy GlobalAdaptive has increased more than \
\\!\\(\\*RowBox[{\\\"2000\\\"}]\\) times. The global error is expected to \
decrease monotonically after a number of integrand evaluations. Suspect one \
of the following: the working precision is insufficient for the specified \
precision goal; the integrand is highly oscillatory or it is not a \
(piecewise) smooth function; or the true value of the integral is 0. \
Increasing the value of the GlobalAdaptive option MaxErrorIncreases might \
lead to a convergent numerical integration. NIntegrate obtained \
\\!\\(\\*RowBox[{\\\"3.375773951740421`*^-8\\\"}]\\) and \
\\!\\(\\*RowBox[{\\\"3.657885544349867`*^-13\\\"}]\\) for the integral and \
error estimates.\"",2,6,2,21519474957897888619,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783265912224*^9}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","slwcon",
   "\"Numerical integration converging too slowly; suspect one of the \
following: singularity, value of the integration is 0, highly oscillatory \
integrand, or WorkingPrecision too small.\"",2,6,3,21519474957897888619,
   "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.8397832660913877`*^9}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","eincr",
   "\"The global error of the strategy GlobalAdaptive has increased more than \
\\!\\(\\*RowBox[{\\\"2000\\\"}]\\) times. The global error is expected to \
decrease monotonically after a number of integrand evaluations. Suspect one \
of the following: the working precision is insufficient for the specified \
precision goal; the integrand is highly oscillatory or it is not a \
(piecewise) smooth function; or the true value of the integral is 0. \
Increasing the value of the GlobalAdaptive option MaxErrorIncreases might \
lead to a convergent numerical integration. NIntegrate obtained \
\\!\\(\\*RowBox[{\\\"1.924972439256785`*^-8\\\"}]\\) and \
\\!\\(\\*RowBox[{\\\"1.8698032241020584`*^-13\\\"}]\\) for the integral and \
error estimates.\"",2,6,4,21519474957897888619,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783274332556*^9}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","slwcon",
   "\"Numerical integration converging too slowly; suspect one of the \
following: singularity, value of the integration is 0, highly oscillatory \
integrand, or WorkingPrecision too small.\"",2,6,5,21519474957897888619,
   "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783275296496*^9}],

Cell[BoxData[
 TemplateBox[{
  "General","stop",
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"slwcon\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"",2,6,6,21519474957897888619,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783275315724*^9}],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","eincr",
   "\"The global error of the strategy GlobalAdaptive has increased more than \
\\!\\(\\*RowBox[{\\\"2000\\\"}]\\) times. The global error is expected to \
decrease monotonically after a number of integrand evaluations. Suspect one \
of the following: the working precision is insufficient for the specified \
precision goal; the integrand is highly oscillatory or it is not a \
(piecewise) smooth function; or the true value of the integral is 0. \
Increasing the value of the GlobalAdaptive option MaxErrorIncreases might \
lead to a convergent numerical integration. NIntegrate obtained \
\\!\\(\\*RowBox[{\\\"3.375773951740421`*^-8\\\"}]\\) and \
\\!\\(\\*RowBox[{\\\"3.657885544349867`*^-13\\\"}]\\) for the integral and \
error estimates.\"",2,6,8,21519474957897888619,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.839783281839622*^9}],

Cell[BoxData[
 TemplateBox[{
  "General","stop",
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"eincr\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"",2,6,9,21519474957897888619,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.819749628966176*^9, 3.819754750088048*^9, 
  3.819941622240633*^9, 3.820036739233984*^9, 3.820073743895574*^9, 
  3.820462278501912*^9, 3.822269604833716*^9, 3.8222730398749733`*^9, 
  3.8397715723781853`*^9, 3.83978328185857*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.`", ",", "0.42976856069869773`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.5`", ",", "0.4301728831764646`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1.`", ",", "0.4313894925968736`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1.5`", ",", "0.43343587570984693`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2.`", ",", "0.43634521772457097`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2.5`", ",", "0.440161220316608`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3.`", ",", "0.44494485892911256`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3.5`", ",", "0.4507742502406256`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4.`", ",", "0.45775042547975753`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4.5`", ",", "0.46599930469488`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"5.`", ",", "0.47567333180895166`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"5.5`", ",", "0.48696626225306827`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6.`", ",", "0.5001191657633984`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6.5`", ",", "0.5154239973062508`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7.`", ",", "0.5332385246847466`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7.5`", ",", "0.5540034321912195`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8.`", ",", "0.5782421753376443`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8.5`", ",", "0.6065524969089802`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"9.`", ",", "0.6395690862387072`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"9.5`", ",", "0.6778449381804675`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"10.`", ",", "0.7215939796756664`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"10.5`", ",", "0.7702198188044326`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11.`", ",", "0.8216618046667364`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11.5`", ",", "0.8719775472042627`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"12.`", ",", "0.9160577559827239`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"12.5`", ",", "0.9498554258680351`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"13.`", ",", "0.9724582616233911`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"13.5`", ",", "0.9858705186846288`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"14.`", ",", "0.9931246131477809`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"14.5`", ",", "0.9967909316669397`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"15.`", ",", "0.9985533374127638`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.819750361046006*^9, 3.819755564066687*^9, 
  3.819942366853595*^9, 3.8200375493329077`*^9, 3.8200744442148314`*^9, 
  3.820462980264586*^9, 3.822270313654598*^9, 3.822273798443659*^9, 
  3.83977230444436*^9, 3.83978396481493*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"calcsuppfactors", "=", 
  RowBox[{"{", 
   RowBox[{
   "0.4313894925968736", ",", "0.43634521772457097", ",", 
    "0.44494485892911256", ",", "0.45775042547975753", ",", 
    "0.47567333180895166", ",", "0.5001191657633984", ",", 
    "0.5332385246847466", ",", "0.5782421753376443", ",", 
    "0.6395690862387072", ",", "0.7215939796756664", ",", 
    "0.8216618046667364", ",", "0.9160577559827239", ",", 
    "0.9724582616233911"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.839774275981242*^9, 3.8397743447265577`*^9}, {
  3.8397743820160646`*^9, 3.839774521720811*^9}, {3.8397745597752542`*^9, 
  3.839774568805504*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.4313894925968736`", ",", "0.43634521772457097`", ",", 
   "0.44494485892911256`", ",", "0.45775042547975753`", ",", 
   "0.47567333180895166`", ",", "0.5001191657633984`", ",", 
   "0.5332385246847466`", ",", "0.5782421753376443`", ",", 
   "0.6395690862387072`", ",", "0.7215939796756664`", ",", 
   "0.8216618046667364`", ",", "0.9160577559827239`", ",", 
   "0.9724582616233911`"}], "}"}]], "Output",
 CellChangeTimes->{3.839774576022743*^9, 3.839784411436907*^9}]
}, Open  ]],

Cell["data :", "Text",
 CellChangeTimes->{{3.8397906014304256`*^9, 3.839790616475357*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"densitymean", "=", 
  RowBox[{
   RowBox[{"0.63", "/", "0.85"}], "*", 
   RowBox[{"{", 
    RowBox[{
    "2.0007", ",", "1.9238", ",", "1.8166", ",", "1.6459", ",", "1.4441", ",",
      "1.2487", ",", "1.0179", ",", "0.7942", ",", "0.5888", ",", "0.3958", 
     ",", "0.2568", ",", "0.1636", ",", "0.0944"}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.8397739186590033`*^9, 3.8397739186794643`*^9}, {
  3.839773967374426*^9, 3.839773979723648*^9}, {3.839774010656438*^9, 
  3.839774054829607*^9}, {3.839784277035159*^9, 3.83978429076626*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1.4828717647058824`", ",", "1.4258752941176471`", ",", 
   "1.3464211764705882`", ",", "1.2199023529411765`", ",", 
   "1.0703329411764706`", ",", "0.9255070588235293`", ",", 
   "0.7544435294117647`", ",", "0.5886423529411765`", ",", 
   "0.436404705882353`", ",", "0.29335764705882356`", ",", 
   "0.19033411764705882`", ",", "0.1212564705882353`", ",", 
   "0.06996705882352941`"}], "}"}]], "Output",
 CellChangeTimes->{3.8397842919657917`*^9, 3.839784415574245*^9}]
}, Open  ]],

Cell["data :", "Text",
 CellChangeTimes->{{3.8397906014304256`*^9, 3.839790616475357*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"pbmean", "=", 
  RowBox[{
   RowBox[{"0.63", "/", "0.85"}], "*", 
   RowBox[{"{", 
    RowBox[{
    "1.0216", ",", "1.0483", ",", "0.9018", ",", "0.7639", ",", "0.6676", ",",
      "0.6264", ",", "0.4879", ",", "0.4316", ",", "0.3215", ",", "0.2034", 
     ",", "0.1462", ",", "0.1288", ",", "0.0581"}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.839779246045546*^9, 3.8397792614871273`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.7571858823529413`", ",", "0.7769752941176471`", ",", 
   "0.6683929411764706`", ",", "0.566184705882353`", ",", 
   "0.4948094117647059`", ",", "0.46427294117647055`", ",", "0.36162`", ",", 
   "0.3198917647058824`", ",", "0.23828823529411766`", ",", 
   "0.15075529411764707`", ",", "0.10836`", ",", "0.09546352941176471`", ",", 
   "0.04306235294117647`"}], "}"}]], "Output",
 CellChangeTimes->{3.839779266368516*^9, 3.839784551158935*^9}]
}, Open  ]],

Cell["theory :", "Text",
 CellChangeTimes->{{3.839790633664722*^9, 3.8397906419220867`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"th", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"densitymean", "[", 
      RowBox[{"[", "n", "]"}], "]"}], " ", 
     RowBox[{"calcsuppfactors", "[", 
      RowBox[{"[", "n", "]"}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "13"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8397748646237297`*^9, 3.839774905089356*^9}, {
  3.8397788502754602`*^9, 3.839778853571865*^9}, {3.839784503697523*^9, 
  3.83978450590372*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.6396952981627011`", ",", "0.6221738656598514`", ",", 
   "0.5990831804238756`", ",", "0.5584108211025809`", ",", 
   "0.5091288362742864`", ",", "0.46286381816696004`", ",", 
   "0.40229835458148266`", ",", "0.3403778346605753`", ",", 
   "0.2791109589714483`", ",", "0.21168511200946605`", ",", 
   "0.15639027459553329`", ",", "0.11107793034544398`", ",", 
   "0.06804004439443095`"}], "}"}]], "Output",
 CellChangeTimes->{3.8397749162886887`*^9, 3.839778856323246*^9, 
  3.839784514113583*^9, 3.839784555397292*^9}]
}, Open  ]]
},
WindowSize->{808, 755},
WindowMargins->{{12, Automatic}, {Automatic, 24}},
FrontEndVersion->"11.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 28, \
2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 635, 9, 30, "Text"],
Cell[1196, 31, 854, 12, 30, "Text"],
Cell[2053, 45, 2165, 54, 138, "Input"],
Cell[4221, 101, 541, 7, 30, "Text"],
Cell[4765, 110, 533, 17, 32, "Input"],
Cell[CellGroupData[{
Cell[5323, 131, 428, 6, 32, "Input"],
Cell[5754, 139, 525, 8, 32, "Output"]
}, Open  ]],
Cell[6294, 150, 97, 1, 30, "Text"],
Cell[6394, 153, 4059, 104, 285, "Input"],
Cell[10456, 259, 200, 4, 30, "Text"],
Cell[10659, 265, 1911, 54, 117, "Input"],
Cell[CellGroupData[{
Cell[12595, 323, 450, 12, 32, "Input"],
Cell[13048, 337, 581, 11, 42, "Message"],
Cell[13632, 350, 1131, 18, 113, "Message"],
Cell[14766, 370, 583, 11, 42, "Message"],
Cell[15352, 383, 1132, 18, 113, "Message"],
Cell[16487, 403, 581, 11, 42, "Message"],
Cell[17071, 416, 554, 10, 24, "Message"],
Cell[17628, 428, 1131, 18, 113, "Message"],
Cell[18762, 448, 552, 10, 24, "Message"],
Cell[19317, 460, 2783, 68, 159, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22137, 533, 649, 13, 96, "Input"],
Cell[22789, 548, 512, 10, 54, "Output"]
}, Open  ]],
Cell[23316, 561, 90, 1, 30, "Text"],
Cell[CellGroupData[{
Cell[23431, 566, 566, 11, 75, "Input"],
Cell[24000, 579, 512, 10, 54, "Output"]
}, Open  ]],
Cell[24527, 592, 90, 1, 30, "Text"],
Cell[CellGroupData[{
Cell[24642, 597, 416, 9, 54, "Input"],
Cell[25061, 608, 486, 9, 54, "Output"]
}, Open  ]],
Cell[25562, 620, 92, 1, 30, "Text"],
Cell[CellGroupData[{
Cell[25679, 625, 500, 13, 32, "Input"],
Cell[26182, 640, 563, 11, 54, "Output"]
}, Open  ]]
}
]
*)

